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I. INTRODUCTION 

Many years ago, Steven Weinberg mentioned to me 
that he was disconcerted that none of the classic text- 
books on superconductivity would explain the phe- 
nomenon in terms of the Higgs mechanism jlj] for the 
electromagnetic gauge field. This concern is of course 
very well justified, and it was most likely with this con- 
cern in mind that Weinberg has included a section on 
superconductivity in his treatment of spontaneous sym- 
metry breaking and the Higgs mechanism in his series 
of volumes entitled The quantum theory of fields 0- 
When I was asked recently to present some lectures on 
superconductivity, I opened his book expecting to find 
a particularly lucid exposition of this in condensed mat- 
ter physics rarely emphasized perspective. I found the 
exposition I was looking for, but to my surprise, build 
around the following statement: A superconductor is sim- 
ply a material in which electromagnetic gauge invariance 
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is spontaneously broken. What Weinberg means with 
this statement is just that the electromagnetic gauge field 
"acquires a mass" due to the Higgs mechanism in a su- 
perconductor, as particle physicists often speak of spon- 
taneously broken gauge invariance interchangeably with 
the Higgs mechanism. 

Nonetheless, I am not perfectly at ease with the above 
statement, which is, by the way, by no means specific to 
Weinberg's exposition, but widely believed and accepted. 
While it is obvious that Weinberg fully understands the 
matter, the statement may still be misleading to a young 
student who is learning the subject for the first time. 
The problem is that the statement is, if one takes it liter- 
ally, not correct: gauge invariance cannot spontaneously 
break down as a matter of principle, and in particular is 
not broken in a superconductor, as I will explain in the 
following section. 

This paper is organized as follows. In Section m we 
discuss the statement quoted above including the dan- 
ger which may result from a literal interpretation of it in 
depth. In particular, we show that the BCS ground state 
is, in contrast to statements made in the literature, fully 
gauge invariant. The crucial ingredient often omitted is 
that gauge transformations involve, in addition to the 
standard transformation of gauge fields, local phase ro- 
tations of both creation (and annihilation) operators and 
wave functions. In Section IlIII we discuss the nature of 
the order parameter in superconductors, with a partic- 
ular emphasis on finite systems, which always possess a 
unique ground state. The arising subtlety is explained by 
drawing an analogy to quantum antiferromagnets, which 
also possess a unique and rotationally invariant ground 
state for finite systems. In Section llVI we motivate and 
elaborate the effective field theory of a superconductor 
at low temperatures, which contains the theory of a neu- 
tral superfluid as the special case where the charge is set 
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to zero. In particular, we obtain the particle density and 
current as well as the energy and momentum density from 
the physical symmetries of the theory, invariance under 
global U(l) phase rotations of the order parameter and 
invariance under translations in time and space. The 
quest for a consistent definition of the superfluid velocity 
yields a relation between current and momentum densi- 
ties in the superfluid, which in turn requires corrections 
to the effective Lagrangian. Since the density of the su- 
perfluid is essentially the "momentum conjugate" of the 
order parameter phase, Hamilton's equations yield phys- 
ical information not contained in the Euler-Lagrange 
equations; specifically, we obtain a gauge invariant gener- 
alization of Anderson's characteristic equations of super- 
fluidity to the case of superconductors. We conclude this 
Section with a brief justification of the classical nature 
of the effective field theory. In Section we discuss the 
phenomenology of superconductors as compared to neu- 
tral superfluids, or, in general terms, the Higgs mecha- 
nism. To begin with, we briefly address the phenomenol- 
ogy of neutral superfluids including vortex quantization, 
and give a general introduction to the Higgs mechanism 
in held theories. We then turn to the phenomenology of 
simply connected superconductors, solve the equations 
of motion, obtain electric and magnetic screening, Lon- 
don's equation, the Bernoulli Hall effect, and the balance 
of the Lorentz force. We demonstrate that the Higgs 
mechanism never corresponds to a spontaneous violation 
of a gauge symmetry, and that it is incorrect to inter- 
pret it in terms of "a mass acquired by the electromag- 
netic gauge field", as the massive field is no longer a 
gauge field. Specifically, we show that the massive vector 
field, which may alternatively be used to describe a (sim- 
ply connected) superconductor, is correctly interpreted 
as a four- vector formed by the chemical potential and the 
three components of the superfluid velocity. We conclude 
this Section with a discussion of the subtle difference be- 
tween the physical invariance of the theory under global 
U(l) rotations of the order parameter phase and gauge 
invariance, which is nothing but a local invariance of our 
description of the system. In the last Section, we review 
a family of "quantum effects" : the quantization of mag- 
netic flux in superconductors, phase slippage, and the 
Josephson effect in both neutral superffuids and super- 
conductors. In these effects, Planck's constant manifests 
itself in the phenomenology through the compactness of 
the order parameter phase field; these effects require ei- 
ther a non-trivial topology or more than one superfluid. 
We derive them from the effective field theory introduced 
in section lTVI and thereby demonstrate that the very few 
assumptions made in motivating the effective theory are 
sufficient to account for them. 



II. GAUGE INVARIANCE 

To illustrate how dangerous the statement quoted in 
the introduction is in the case of superconductivity. 



where we do not only have a description in terms of an 
effective field theory but also a microscopic description 
in terms of model Hamiltonians and trial wave functions, 
I will at first assume the statement was true and take 
it literally. I will pretend to be a student who has just 
learned that electromagnetic gauge invariance is sponta- 
neously violated in a superconductor. Well, what does 
this mean? A spontaneously broken symmetry means 
that the Hamiltonian of a given system in the thermody- 
namic limit is invariant under a given symmetry trans- 
formation [i.e., commutes with the generator(s) of this 
symmetry) while the ground state is not invariant. There 
are many ground states, which transform into each other 
under the symmetry transformations. A classic example 
is ferromagnetism: The Hamiltonian is rotationally in- 
variant, while any particular ground state, specified by 
the direction the magnetization vector points to, is not. 
So if gauge invariance is broken in a superconductor, this 
must mean that the ground state of the superconductor 
does not share the gauge invariance of the Hamiltonian. 
Indeed, a glance at the BCS wave function 

k 

where the coefficients Uk and Vk are chosen real and (/> 
is an arbitrary phase, appears to confirm this picture. 
There are many different ground state wave functions, 
labeled by 4>, which transform into each other under an 
electromagnetic gauge transformation given by 



which is tantamount to taking 

. ^ 2e , 
0^(^-f — A. 
Tic 

For simplicity, we have chosen A independent of space- 
time. The electron charge throughout this article is — e. 

Next, I the student ask myself whether these many 
BCS wave functions for different parameters (j) corre- 
spond to physically different states. I know that gauge 
transformations are not physical transformations: gauge 
invariance is an invariance of a description of a system, 
while other symmetries correspond to invariances under 
physical transformations, like rotations or translations, 
which affect the physical state in question. For example, 
if the Hamiltonian for given system (like a ferromagnet) 
is invariant under rotations in space, this implies that 
if we rotate a given eigenstate, we will obtain another 
eigenstate. Depending on whether the original state is ro- 
tationally invariant or not, it will transform into itself or 
into a physically different state. A gauge transformation, 
by contrast, will only transform our description of a sys- 
tem from one gauge to another, without ever having any 
effect on the physical state of the system. Gauge trans- 
formations are comparable to rotations or translations of 
the coordinate system we use to describe a system. An- 
other way of seeing the difference is by noting that it is 
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possible to rotate or translate a superconductor in the 
laboratory, but as a matter of principle not possible to 
gauge transform it. Returning to the superconductor, I 
the student conclude that if the many different ground 
states only differ by a gauge transformation, they cannot 
be physically different. The ground state of a supercon- 
ductor must hence be physically unique. 

In fact, there is another way of looking at the prob- 
lem which appears to confirm this conclusion. A BCS 
superconductor can not only be described in the grand- 
canonical ensemble, where the chemical potential rather 
than the number of particles is fixed, but also in the 
canonical ensemble, where the number of particles or 
pairs is fixed. Following Anderson 0, we can project 
out a (not normalized) state with N pairs from via 

r"d</)e-*^*|V0) (3) 
and obtain (see appendix) 

|V'JV)= Jd^Xi. ..(fx2N'P{Xi -X2) ■■■■■ (p{x2N-l- X2n) ' 

• V| (a;i)?A|(a;2) • • ■ i'\ix2N-i)''J^l{x2N) I ), (4) 

where the real-space creation operator fields 'iptix) are 
simply the Fourier transforms of the momentum-space 
creation operators c^^, 

i^Ux) = -L^e-'=-cL, cl = Jd^xe^^-^lix). 

, _ , (5) 

The wave function for each of the individual pairs, which 
only depends on the relative coordinate, is (up to a nor- 
malization) given by 

This form nicely illustrates that all the pairs have con- 
densed into the same state, which is the essence of su- 
perfluidity. As (p{x) is uniquely determined for a given 
Hamiltonian, the ground state of the superconductor Q 
once more appears to be unique and non-degenerate. 

So far the students train of thought. The conclu- 
sion reached is of course completely wrong: a superfluid, 
and in particular a superconductor, is characterized by 
a spontaneously broken symmetry, and, at least in the 
thermodynamic limit, there are many degenerate ground 
states. There are several mistakes in the students analy- 
sis. The first is his literal interpretation of the statement 
quoted in the inroduction. In fact, a gauge symmetry 
cannot spontaneously break down as a matter of princi- 
ple, since it is not a physical symmetry of the system 
to begin with, but merely an invariance of description 
The only way to violate a gauge symmetry is by choosing 
a gauge, which again has only an effect on our descrip- 
tion, but not on the physical system itself. 



In particular, the BCS ground state does not violate 
gauge invariance, even though statements to the contrary 
have been made in the literature. The apparent contra- 
diction with and Q can be resolved by recalling that 
a gauge transformation only affects our description of the 
system, and is analogous to a rotation of the coordinate 
system we use in the example of a ferromagnet: if we 
rotate the coordinate system accordingly, a ground state 
with the magnetization vector pointing in the ^-direction 
in the original coordinate system will "transform" into a 
state with the magnetization vector pointing in the x- 
direction in the new coordinate system, while the physi- 
cal state has not been affected at all. So while the BCS 
wave function may look different in a different gauge, the 
state itself will remain the same. 

It is worthwhile to rephrase this statement in equa- 
tions. To begin with, let us consider a (relativistic quan- 
tum) field theory. Electromagnetic gauge invariance is 
the invariance of a given theory under U(l) rotations of 
the complex scalar fields which carry the charge: 

V'^(x) — > e*'ic^(^-'?/;1'(x), "ij^ix) — > e~*?»c^'-^^^(x), (7) 

where x denotes spacetime. We use the conventions 
(x^) = {x°,x'^,x'^,x^) = {ct,x,y,z), {Xf,) = g^i/x'', 
1 = 500 = -ffii = -.922 = -533- If the theory con- 
tains gradient terms in these fields (as it usually does), 
gauge invariance demands that they are minimally cou- 
pled to a U(l) gauge field, i.e., the gradient terms must 
enter the Lagrangian as 

(d^ + i-^Af,{x)'jij\x) or (^d^, ~ i-^Af,{x)jTl>{x), 

where (9^) = (d/dx^) — (i9t,V). The gauge field 
(Af^) = (<I>, —A) must transform according to 

A,,{x)^A,,{x)-df,A{x). (8) 

The statement that the theory is gauge invariant simply 
means that the Lagrangian is invariant under the com- 
bined transformation ||7J) and ©. It is not a physical 
invariance, but an invariance of description, as it only 
amounts to a reparamctrization of fields. 

The concept of gauge invariance is implemented in a 
very similar way in non-relativistic quantum mechanics, 
where the gauge field A is no longer considered a dy- 
namical variable, but an externally applied vector po- 
tential, and we usually do not describe a system by a 
Lagrange density, but by a Hamiltonian operator and its 
eigenstates. For pedagogical reasons, let us first assume 
a formulation in second quantization. Electromagnetic 
gauge invariance means once again that the description 
is invariant under U(l) rotations of the particle creation 
and annihilation operator fields , 

^lix) ^ e*i''("Vi(a^), Mx) ^ e~'fc^^'''>Mx). 

.... (9) 
The kinetic part of the Hamiltonian will again contain 
gradient terms in the operator fields, which once again 
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must be minimally coupled to the electromagnetic gauge 
field. For example, the standard kinetic Hamiltonian for 
a quadratic dispersion 

(10) 

is obviously invariant under ^ provided we transform 
the gauge field simultaneously according to 

A{x) ^ A{x) +VK{x). (11) 

Let us now turn to the gauge transformation properties of 
the eigenstates. Consider a general N electron eigenstate 

|(/7)= ^ j (f'Xi...SxNip{xi...XN](Ti...aN)- 
(Ti...aN 

^lix,)...^ljx^)\0). (12) 

The state is invariant under provided we transform 
the wave function according to 

ip{xi . . .X]^;ai . . .(7n) -> 

N 

Y[ e-'7fe'^(^^' ifixi ...XN\ai... gn). (13) 

This already illustrates the statement phrased in words 
above: A gauge transformation leaves physical states in- 
variant. This is just not obvious in every formulation. 
If we formulate a problem in non-relativistic quantum 
mechanics in first quantization, a gauge transformation 
will only amount to (|ll|l and H13|l . as we do not even 
introduce the operator fields -0^ and As A{x) im- 
plements an externally applied magnetic field, we must 
choose a gauge in order to obtain explicit expressions for 
the Hamiltonian and the eigenstates. The vector poten- 
tial and the wave functions will have different functional 
forms in different gauges. The gauge rotations © of the 
particle creation and annihilation operators, by contrast, 
only amount to a local change of variables; we could write 

^t(a;)^Vt'(a^) = e'^^^(-)Vi(a;), 
^„{x)^^J{x) = e-'tc^^'^^ ^4x), 

and then simply omit the primes. This part of the gauge 
transformation is often omitted as a choice of convention. 

In the case of a BCS superconductor, such a convention 
would be all but propitious, as it would suggests that the 
ground state is not gauge invariant. The apparent con- 
tradiction in the students train of thought is immediately 
resolved as one uses the full and correct prescription for 
a gauge transformation, 

cL^e^^c^cl, cb^cb-^^A, (14) 

where the transformation of the phase cf) is the equivalent 
of Ijl^ll above. Then the BCS ground state 

[](ufe + «fce^^4^cL^^) |0), 

k 



is evidently gauge invariant; it is merely the label in 
\ip^) which will be adjusted under a gauge transforma- 
tion. The transformation (f> (p — ||A is also required 
for the classical (or Ginzburg-Landau) order parameter 
field '5* (a;), which is given by the expectation value of 
the operator field 

i>Hx) = i:\ix) ^Plix), (15) 

to have to the correct gauge transformation properties. 
For the BCS ground state, 

^*{x) = {^4 ^P\{x) ^lix) \^^) = 1 ^ vlu^e-^''' (16) 

k 

will transform as a field of charge — 2e under H14|l . 

**(a;) -> e'ii^(^)**(a;). (17) 

This is the physically correct prescription. When we cou- 
ple 'i'*{x) minimally to the electromagnetic gauge field, 
as required by ((T7|l . we obtain the correct effective field 
theory description of superconductivity. This theory dis- 
plays the Higgs mechanism and yields London's equation. 
(By contrast, if we were to adhere to (O, 5** (a;) would 
be invariant, could not be coupled to the electromagnetic 
gauge field, and no sensible effective field theory could be 
formulated.) 

III. ORDER PARAMETER CONSIDERATIONS 

Before proceeding further with the Higgs mechanism, I 
would like to return to the students train of thought and 
explain what is wrong with his conclusion drawn from 
the BCS wave function in position space. The problem 
here is that while the ground state is indeed unique for 
a finite system, there are many degenerate states in the 
thermodynamic limit, which correspond to different num- 
bers of particles. To understand this issue in depth, it 
is best to first recall how rotational symmetry is sponta- 
neously violated in ferromagnets and antiferromagnets. 
As a minimal model, we consider a three dimensional 
cubic lattice of spins with spin quantum number S and 
assume the Heisenberg Hamiltonian |lO|| 

HJ = JY,S^SJ (18) 

where the sum extends over all nearest-neighbor bonds 
(«, j) and J < (J > 0) for a ferromagnet (an antifcrro- 
magnet). 

In the case of a ferromagnet, the order parameter is 
given by the total spin operator 

5tot = ^5^, (19) 

i 

where the sum extends over all lattice sites. It commutes 
with the Hamiltonian, 

[Hj,Stot]^0, (20) 
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and it is hence possible to choose simultaneous eigen- 
states of the Hamiltonian and the order parameter. In 
other words, the degenerate eigenstates of the order pa- 
rameter corresponding to all possible directions the mag- 
netization vector 

M = (Stot) 

can point to, are simultaneously degenerate eigenstates 
of the Hamiltonian. (For a ferromagnet with N spins, all 
the spins align and the ground states are just the states 
with maximal total spin Stot = NS.) The ground state of 
the Hamiltonian is vastly {i.e., 2S'tot + l fold) degenerate 
even if the system is finite. 

The situation is different in the case of the antiferro- 
magnet. The order parameter is given by the Neel vector, 
which in operator form is given by 



(21) 



where A and B denote the two sublattices of the (bipar- 
tite) cubic lattice. It does not commute with the Hamil- 
tonian: 



^0. 



This implies that we cannot choose simultaneous eigen- 
states for the Hamiltonian and the order parameter. In 
fact, a theorem due to Marshall [ij states that the 
ground state for even is unique and a spin singlet, 
or in other words, rotationally invariant. (It is possible 
to choose simultaneous eigenstates of Hj and S'tot as (|20|) 
holds independently of the sign of J.) The classical Neel 
order parameter, 

N = (N) 

will vanish for any finite system. This is not to say that 
there is no order for a finite system; it just manifests it- 
self only through long-range correlations in the staggered 
spin-spin correlation function: 



± const. 



-J 



where the -I- sign applies for i and j on the same sub- 
lattice, the — sign for i and j on different sublattices, 
and i~j oo is understood to denote a very large sep- 
aration within (the finite volume of) the system. As we 
approach the thermodynamic limit, the difference in en- 
ergy between the lowest singlet and lowest eigenstates for 
•S'tot = 1,2,3... vanishes, and the ground state becomes 
degenerate (see Fig. ^). These degenerate states can 
now be classified by the directions of the Neel vector N, 
and the spontaneous breakdown of rotation symmetry is 
evident. 

The situation in superconductors is analogous to the 
antiferromagnet: The (operator valued) order parameter 
(|15|l does not commute with the BCS Hamiltonian for 
any finite system even if we work in the grand-canonical 
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FIG. 1: In antiferromagnets (a) and superconductors (b), the 
ground state is unique for finite systems but degenerate in the 
thermodynamic limit. 



ensemble, and the ground state for any finite volume will 
have a well defined particle number. The difference in 
energy between a system with NoTNztlorNzL2 etc. 
pairs, however, will vanish in the thermodynamic limit 
(see Fig.^), and the many degenerate ground states 
can be classified by the phase cj) of the (classical) order 
parameter 



(22) 



The broken symmetry is of course also present in a system 
with a fixed number of particles, but like in the case of 
the antiferromagnet, only as a long-range correlation of 
the (operator valued) order parameter field: 



{^^{x)'^{y)) -> const. 



y| 



(23) 



where ^(cc) — '>pi{x)ip-^(x) is simply the hermitian conju- 
gate of ^^{x). This correlation is referred to as off diago- 
nal long-range order (ODLRO) 0|. This type of order is 
characteristic to all superfluids, whether charged (like a 
superconductor) or neutral (like liquid helium), whether 
fermionic (like a superconductor or '^He) or bosonic (like 
^He). For a bosonic superfluid, the (operator valued) or- 
der parameter ^^{x) and its hermitian conjugate ^{x) 
no longer create or annihilate a pair of fermions, but sim- 
ply create or annihilate a single boson (like a ^He atom) . 

The ODLRO is already evident from the position space 
wave function Q : Since all the pairs have condensed into 
the same quantum state, which is translationally invari- 
ant as it does not depend on the center-of-mass coordi- 
nates of the pairs, we expect to obtain a finite overlap 
with the original ground state if we rather clumsily (i.e., 
via '^{y)) remove a pair of particles at some location y 
and equally clumsily {i.e., via ^'''(a;)) recreate it at a dis- 
tant location a;. In a superfluid or superconductor with a 
fixed number of particles, the phase will align over the 
entire system, like the direction of the staggered magne- 
tization or Neel vector will align in an antiferromagnet. 
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To illustrate the significance of the phase once more, 
let us consider a large (but finite) superconductor A, and 
describe it as a combination of two superconductors B 
and C: 





If we label the ground states of each superconductor by 
its phase, we can obviously write 

as the phase </> of the order parameter will align over the 
entire system. If we now transform to a description in 
terms of fixed numbers of pairs Na for each supercon- 
ductor. 



d(j) e 



'iNa4> 



where a can be A, B, or C, the ground state of A is no 
longer a direct product of the ground states of B and C: 

no matter how we choose A^b, as the phases no longer 
align. So while it is possible to describe a superconductor 
in a canonical ensemble (i.e., with fixed particle number), 
it is highly awkward to do so. It is comparable to a 
description of an antiferromagnet with long-range Neel 
order in terms of an overall singlet ground state of the 
system. 

The most significant difference between an antiferro- 
magnet and a superfluid or superconductor with regard 
to the order parameter is that the broken rotational sym- 
metry in the former case is much more evident to us, as 
all the macroscopic objects in our daily life experience 
violate rotational symmetry at one level or another. In 
particular, the structure of the material in which anti- 
ferromagnetic order occurs provides us already with a 
reference frame for the direction the Neel order parame- 
ter may point to. In the case of the superconductor, we 
need a second superconductor to have a reference direc- 
tion for the phase, and an interaction between the order 
parameter in both superconductors to detect a relative 
difference in the phases. (In practice, such an interaction 
may be accomplished by a pair tunneling or so-called 
Josephson junction.) The interference experiments will 
of course only be sensitive to the relative phase, and not 
the absolute phase in any of the superconductors, as all 
phases can, as a matter of principle, only be specified 
relative to some reference phase. In principle, the same 
is true for rotational invariance, but in this case the fixed 
stars provide us with a reference frame we perceive as 
"absolute" . 

We may conclude at this point that in a superfluid or 
superconductor, a symmetry is spontaneously violated. 



but this symmetry is not gauge invariance, but global 
U(l) phase rotation symmetry. This is already evident 
from the fact that the discussion above made no reference 
to whether the order parameter field '^'^{x) is charged or 
not, and equally well applies to neutral superfluids, where 
^'^{x) carries no charge. 

There is, however, a very important difference between 
these two cases. If the order parameter field is neutral, 
the excitation spectrum of the system contains a gapless 
(or in the language of particlephysics "massless" ) mode, 
a so-called Goldstone boson [y, which physically corre- 
sponds to very slow spatial variations in the direction (as 
for the case of broken rotational invariance) or phase (as 
for the case of a superfluid) of the classical order parame- 
ter fleld. If the order parameter fleld is charged, however, 
it couples to the electromagnetic gauge field, and the 
Goldstone boson is absent due to the Higgs mechanism. 
The physical principle underlying the mechanism was dis- 
covered by Anderson 13j in the context of superconduc- 
tivity: as the electromagnetic interaction is long-ranged, 
the mode corresponding to very slow spatial variations 
in the phase 4> of the superconducting order parameter, 
which implies currents by the equation of motion and 
hence also variations in the density of the superfluid by 
the continuity equation, acquires a gap (or "mass" ) given 
by the plasma frequency. 



IV. EFFECTIVE FIELD THEORY 

Most of the phenomenology of superfluidity or super- 
conductivity can be derived from a simple effective field 
theory, which in the latter case displays the Higgs mech- 
anism. It is probably best to turn directly to the low- 
energy effective Lagrangian for the superconductor, as 
it contains the superfluid as the special case where the 
coupling e* of the order parameter to the electromagnetic 
gauge fleld is set to zero. To motivate the Lagrangian, re- 
call first the Ginzburg-Landau expansion of the free 
energy density in terms of the order parameter (which is 
now normalized differently from (|16|l above) in the vicin- 
ity of the critical temperature Tc, where the transition 
between normal and superconducting phases occurs: 



1 



2m* 



-ihV H A{x))^{x) 



+ a{T)\^{x)[^ 



\b{T)\^>{xt + ^B{xf 

Z OTT 



(24) 



where m* and — e* = — 2e are the effective mass and 
charge of the electron pairs, respectively, and J5 = V x A 
is the magnetic field. The material parameter a{T) 
changes sign to become negative as we pass through the 
transition from above, while h{T) has to remain posi- 
tive. Minimizing the free energy in the superconducting 
phase yields that (i) the gradient term must vanish, (ii) 
|*(a;)P = -a/b, and (in) B = Q. This means that the 
amplitude ^'o of the order parameter "^{x) = ^I'oe*'^ has 
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to be fixed while the phase 0, which labels the many de- 
generate ground state configurations, can be arbitrary as 
long as the variation over the sample is given by 

V(/. = -f A, 
fie 

which implies 4>{x) — const, if we choose the gauge 
A{x) = 0. In the vicinity of the transition, we may 
treat ^' as a small parameter, which implies that the ex- 
pansion (|24l) provides us with a complete description of 
the system at the level of thermodynamics. 

The Ginzburg-Landau expansion is also helpful in mo- 
tivating the low energy effective Lagrange density at low 
temperatures. To begin with, we may assume that since 
the amplitude fluctuations are massive, they do not en- 
ter in the low energy description. Taking |^'(a;)| to be 
constant, the free energy density above reduces to a con- 
stant, an electromagnetic field contribution, and 



fn 



2ra' 



-A{x) 



(25) 



where Tig = 1 5*0 P is a phenomenological parameter which 
depends on the material and the temperature. It has the 
dimension of a density and is equal to the density of the 
superfluid in the absence of currents and inhomogeneities 
at T = 0, as we shall see below. It is usually referred to as 
the superfluid density, but it would be more appropriate 
to use the superfluid stiffness rig /to* as a parameter in- 
stead (2j ■ We will also see below that Galilean invariance 
of the superfluid implies that to* is the bare mass of the 
superfluid particles, i.e., to* — 2toc for Cooper pairs |l5j| . 

We take H25f) to be part of the potential energy in 
the effective Lagrange density for the superfluid. The 
remaining contribution arises from the coupling of the 
charge of the superfluid to the electrostatic potential 
^(x), which is in leading order given by 

/ci = -nse*$(x). 

This term is usually not included in the free energy of the 
superconductor, as it is always canceled off by another 
such term with opposite sign arising from the uniform 
positive background charge. It is essential to our effective 
field theory here, however, as it is part of the Lagrange 
density for the superfluid, while the uniform background 
charge is accounted for by another Lagrange density 



(26) 



where x — {ct, x) denotes spacetime. Note that /d is not 
invariant under (time dependent) gauge transformations 



0(a;) ^ (f>{x) - j-Hx), 

^{x) ^ <^{x) - -dtA{x), 
c 

A{x) A{x)+VA{x). 



(27) 



We now turn to the kinetic energy term in the effective 
Lagrange density. The simplest gauge invariant Lagrange 
density containing both potential energy terms above is 
given by 



(28) 



Csix) = - n^{hdt(t){x) - e*$(x) 

■{h\/<j){x) + —A(x) 
\ c 



2m* 



This Lagrange density, however, cannot be complete. 
The only term containing a time derivative, dt(t>{x), ap- 
pears as a total derivative (here time derivative of 0) 
and hence does not affect the Euler-Lagrange equations 
of motion. It is nonetheless of physical significance, as it 
both ensures gauge invariance and accounts for the lead- 
ing contribution to the particle density, as we will see 
below. 

To obtain a second order time derivative term, recall 
that the characteristic feature of a neutral {i.e., e* = 0) 
superfluid is that the only excitation at low energies is a 
sound wave with a linear dispersion 



uj{k) — v\k\ 



(29) 



where k is the wave number and v is the velocity of sound 
in the fluid. As we wish the effective Lagrange density 
for the superfluid both to be gauge invariant and to yield 
(|29|l as an equation of motion for e* = 0, we arrive at 

Csix) = -n,(^hdtcbix)-e*<i>{x)^ (30) 
+ £^|^(/i5t0(x)-e*ci>(a;))'-(w0(a;)-f^A(x))' 
With 



(31) 



where (9^) = (^9t,V) and (^^) = {<^,-A), the La- 
grange density may also be written 



Cs = -cnsDo<, 



2m* 1 ti^ 



(i?o0)'-(A0)n, (32) 



where i — 1,2, 3. The total Lagrangian of the system is 
given by 



L = jd^x{Cs{x)+Cb{x) + C,n.{x)}, 



(33) 



where 



Ccn. = --r^A-^"" with Pf^- = 9^^- - duA,, (34) 

denotes the standard Maxwell Lagrange density for elec- 
tromagnetism. 

The astonishing feature is now that this simple La- 
grangian for the compact U(l) fleld 4>{x) (compact since 
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the values (p and + 27r describe the same physical state 
and hence must be identified) coupled to the electromag- 
netic gauge field accounts for all the essential features of 
superfluidity or superconductivity. There are also impor- 
tant corrections to it, but we will discover them automat- 
ically as we proceed. 

In order understand the physical content of the La- 
grangian, it is highly instructive to study its symmetries, 
in particular particle number conservation and invariance 
under translations in space and time. We wish our anal- 
ysis to apply both to the case of a neutral and a charged 
superfluid. In the former case, the theory is no longer 
invariant under a local U(l) gauge transformation (as 
the electromagnetic gauge transformation (0 reduces to 
the identity transformation for e* = 0) but still invariant 
under a global U(l) rotation 



or 



4>{x) (j){x) + A, 



(35) 



where A is independent of spacetime. Physically, this 
symmetry corresponds to particle (or Cooper pair) num- 
ber conservation. According to Noether's theorem if 
under a given transformation the Lagrange density only 
changes by a total derivative. 



dX 



d^Ff^ix) 



A=0 



there is a conserved current associated with this symme- 
try: 



J'^{x) = const, 
where 



Sid^q>{x)) 



D(j){x)- F^ix), \ (36) 



D(j){x) 



d4){x, A) 



dX 



A=0 



Current conservation means 9^ — 0. Since H35|l yields 
F'^{x) — and D(j>{x) — 1, the particle four-current 
( J^) = {cp, J) is given by 



5L, 



h5{d^4>{x)) 5{D^(p{x)y 



(37) 



where we have chosen the normalization such that the 
electric current equals the charge times the particle cur- 
rent: 



-e*J''{x). 



The Lagrange density H30|l yields for the particle density 
1 SLs 



p{x) = - 



ns{dt^{x)) 

ris 1 
m 



^(^hdt^{x)-e*^x) 



(38) 
(39) 



and for the particle current 
J(x) = -- 



hS{\'<p{x)) 

^(hV^{x) + -A{x) 
m* V c 



(40) 



(41) 



The corresponding conservation law is just the continuity 
equation 



dtp + \'J = 0. 



(42) 



Note that since our Lagrangian 133|) does not depend 
on (t>{x), but only on derivatives of (j){x), i.e., 



= 0, 



S^x)) 

(|37|l implies that the current conservation law 



-hd^j'^ix) = d, 



5{d^ct>{x)) 



= 



is equivalent the Euler-Lagrange equation for the field 
(j}{x). For a neutral superfluid, we obtain 



(43) 



and hence the dispersion H29() by Fourier transformation. 

The most important implication of l|37|l for the parti- 
cle four-current is, however, that the density p{x) is up 
to a numerical factor equal to the momentum flcld ■n{x) 
conjugate to 0(a;): 



-hp{x) = Tr{x) 



S{dMx)) 



(44) 



We may hence go over to an Hamiltonian formulation, 
and write the Hamiltonian density 



Tisix) = -hp{x)dt(t>{x) - Cs{x), 



(45) 



which is now considered a functional of p{x), di4>{x), 
$(x), and Ai{x), but not dt(j)(x). (In principle. Tig could 
also depend through Cs on 4>{x) and x. Note also that 
(|45|l as the generator of time translations is not invariant 
under time dependent gauge transformations, while the 
equations of motions below are invariant.) The Hamilto- 
nian is of course given by 



Hs^ j d^xHsix). 



(46) 



Hamilton's equations are in analogy to the familiar equa- 
tions 

._ dH{p,q) . ^ dH{p,q) 
dp ' dq 

from classical mechanics given by 

SH, 1 SH^ 



dt(f){x) 



5{-k{x)) h5{p{x)) 



(47) 
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and 



-hdtp{x) = dtTT{x) = di 



5{d,<p{x)) 5{ct>{x)y 



(48) 



With regard to the explicit equations of motion for 
the fields, these equations are equivalent to the Euler- 
Lagrange equation. They provide, however, additional 
information regarding the physical interpretation. To ex- 
tract this information, it is best to study first the other 
conservation laws corresponding to energy and momen- 
tum. 

The theory is invariant under spacetime translations 
X — *■ X — eA, where e is an arbitrary unit vector in space- 
time (e.g. e = (1,0,0,0) or e'^ = Soi^ for a translation 
in time) . The infinitesimal translations are equivalent to 
the field and density transformations ^3 

(j){x) (j){x + eA) = <f){x) + Xe''d^(f>{x), 

Af,{x) Ai,{x + eA) = Af,{x) + Xe^dr^Af^ix), 

C{x) C{x + eA) = C{x) + Xd^ie^Cix)), 

which implies = e'^d^cj), DA^ = e'^d,yA^^ and 

F^(x) = e''£(x). The conserved current associated with 
this symmetry is according to (|36l) given by 



5L 



6L 



e di/A^ £> — €i^Tt[ ^ 



where the canonical energy-momentum tensor T^J^ is 
the sum of the contributions from the superfluid, the 
uniformly charged background, and the electromagnetic 
field: 



where 



rpf-LU 



(49) 



(50) 



The conservation law d^jTlfJ^ = describes energy conser- 
vation for v — Q and momentum conservation for v = i. 
The fi — components of T^J^ correspond to energy and 
momentum densities; in particular, — c Jd^xT^^i gener- 
ates translations in time and J (Px T^^^ translations in 
space. 

In the case of a gauge theory, like the theory of a 
charged superfluid we consider here, it is not possible to 
interpret T^^^^^ as the energy or \T^I^^ as the kinemati- 
cal momentum density of the superfluid. The reason is 
simply that H49|l (and also HSOfl ) is not gauge invariant. 
To circumvent this problem, we simply supplement the 



naive translations by a suitable gauge transformation, 
such that the fields transform covariantly: 

->0 + Ae''(a,0-£^,), 

A^^ A^ + Xe'^id^A^-dM- 

The gauge transformation is hence given by (|27|l with 
A(a;) = Xe'^Ai,{x). This yields the "kinematical" energy 
momentum tensor, with contributions from the super- 
fluid and the electromagnetic field 



(51) 



These expressions are manifestly gauge invariant. Using 
(|37|l with /i = or H38|l . we can write the energy density 
of the superfluid 

T™(a;) = -p{x){hdt<P{x) ~ e*$(x)) - C,{x). (52) 

Note that this is numerically equal to 

T^'^ix) = Hsix) + p{x)e*<^{x). (53) 
Similarly, we can write the momentum density 

-T^\x) = p(x)(w0(x) -I- -A{x)). (54) 
c V c / 

We can use this expression to introduce the superfluid 
velocity Vs{x). In terms of Vs{x), the momentum density 
of the superfluid has to be given by 



-T^\x)=p{x)m*v,{x), 



(55) 



which leads us to define 



m*v^{x) = h\/(j){x) + —A{x) = Di<j){x). (56) 
c 

Since Vs{x) is to be interpreted as a physical velocity, it 
has to transform like a velocity under a Galilean trans- 
formation, 



Vs{x) ~* Vs{x) + U. 



(57) 



The total momentum of the superfluid will hence trans- 
form according to 

i jd^xT^\x) ^ i jd^xT^\x)+ j(fxp{x) m*u, 

(58) 

which implies directly that in a translationally invariant 
system, m* has to be the bare mass of the superfluid 
particles or Cooper pairs . 

It should also be possible to express the particle cur- 
rent in terms of the superfluid velocity. Since the same 
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particles which carry the momentum also carry the cur- 
rent, the particle current has to be given by 



J{x) = p{x)Vsix). 



(59) 



This is almost, but not quite, equivalent to our earlier 
expression (|41ll , as Ug is only the leading contribution to 
p{x). So either H59() with H56() or H41|) is not fully cor- 
rect. To see which one, recall that we have only used 
the general expression for the density 1)38(1 as defined 
through particle number conservation in obtaining l|54|) 
and hence (|56|l and H59() from (|51() . while we have used 
the explicit expression for the Lagrange density ((30(1 in 
obtaining 1(41(1 . In other words, only symmetry consid- 
erations enter in ((54(1 . while ((41(1 depends explicitly on 
the Lagrange density. The expression for the momentum 
density 1(541) . and hence our definition of the superfluid 
velocity 1(561) . is therefore exact, while the expression l(41() 
for the particle current is only an approximation [l9j |. 

The expression for the current, however, will assume 
the exact and physically correct form 1(59(1 if we intro- 
duce suitable corrections to the effective Lagrangian. To 
obtain these, we simply require the Lagrangian to sat- 
isfy 113 



1 



where p{x) and Vs{x) are given by l(M|) and (jSB|l . 

Let us now return to Hamilton's equations, and in par- 
ticular their physical interpretation. With ((53() we may 
rewrite ((471 as 



hdt(l>{x) = - 



5H, 
5{p{x)) 



dTisix) 9T00(a;) 



dp{x) 



dp{x) 



-p{x) +e*<I>(x), 



e*<^{x) 
(66) 



T^\x)^m*.r{x) 



(60) 



where we have used the definition of the chemical po- 
tential. This is one of two equations Anderson J\ refers 
to as "characteristic of superfluidity" . In analogy to the 
definition l(56() of Vs{x), we rewrite it for later purposes 
as 

-p{x) = hdt(j){x) - e*$(a;) = cDo(j){x). (67) 

Taking the gradient and adding ^dtA{x) on both sides 
of ((66() . we obtain 

dt(hV(j){x) + ^Mx)) = -Vpix) - e*E{x) 

= -VAi.,.,„,.(a;), (68) 
where we used the definitions of the electric field, 



1 i ^ SLs 

-U (p = —TO 



c 6{Do(j)) ' 



<5(A0) ' 



(61) 



Upon integration of this equation we find that the La- 
grange density must be of the form 



C.=P[cD, 



1 



■{D. 



2 ' 



(62) 



where P is an arbitrary polynomial. Our superfluid La- 
grange density ((32() will assume this form if we add third 
and fourth order corrections in Dy6\ the full superfluid 
Lagrange density is then given by |20l l2l| 



£s = - n, 



V 2m* / 

-i(ci^o0 + ^(A 



ris 1 



2to* 

It yields for the particle density 
1 <5Ls 



(63) 



p{x) 



c S{Do^) 



1 



Us I / 

Us 7T I CVq^ 

TO* \ 2m 

TO 



and for the particle current 



■(A 



(64) 



J{x) 



<5(A 



p{x)vs{x), 



(65) 



E = -V$ - -dtA, 

c 

and of the electrochemical potential. With 1(56(1 we may 
write 



m* dtVs{x) 



(69) 



The gradient of the electrochemical potential (or chemi- 
cal potential for a neutral superfluid) is usually 01 identi- 
fied with minus the total force on the particles, and II69() 
is referred to as the "acceleration equation". This is, 
however, not quite correct. dtVg in 1(69(1 denotes the time 
derivative in the superfluid velocity field at spacetime x 
(known as "local acceleration" in hydrodynamics), while 
the force on the particles is given by the time deriva- 
tive of the velocity of a given particle in the fluid at x 
("substantial acceleration" in hydrodynamics): 



m* at 



(70) 



Nonetheless, ((69(1 is one of the fundamental equations in 
the phenomenology of superfluidity. It states that if there 
is a gradient in the electrochemical (or chemical) poten- 
tial in a superconductor (or superfiuid), the superfluid 
will be "accelerated" without any frictional damping. On 
the other hand, if the superfluid flow is stationary, the 
electrochemical (or chemical) potential has to be con- 
stant across the superconductor (or superfluid). Since 
a voltmeter measures a difference in the electrochemical 
potential, there cannot be a voltage across a supercon- 
ductor unless the flow is "accelerated" . 
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Let us now turn to Hamilton's second equation (|48|l . 
We first rewrite it as 



-hdtpix) = di 



(71) 



Since the last term in (|53|l . p{x)e*^{x), does not depend 
on the phase 4>ix), we may replace Ti.s{x) in the last term 
in (|71|l by T^'^{x). Integrating the resulting equation over 
the superfluid, discarding a boundary term, and defining 
a "global" derivative with respect to the phase, 

dF[cl,{x)] _ ^.^ F[0(x) + A<^] 



A(j) 



where F[(j){x)] is an arbitrary functional of (pix) and Acf) 
an infinitesimal independent of spacetime, yields 



hdt / (fix p{x) 



hdtN ^ 



d 



d'x '(x) 



(72) 



where N is the number of particles (or pairs) and the 
energy of the superfluid. This is the other "characteristic 
equation" of superfluidity 0- 

This concludes our derivation or motivation of the fun- 
damental equations of superfluidity in the limit of low 
temperatures and low energies. Before turning to the 
phenomenology these equations imply, I would like to di- 
gress briefly and justify one of the implicit assumptions 
made above. The assumption is that we can describe the 
macroscopic quantum phenomena of superfluidity with a 
classical effective field theory, or in other words, that we 
may consider both the phase 4>{x) and its conjugate field 
7r(a;) = —hp{x) as thermodynamic variables. To justify 
this assumption, let us canonically quantize the theory 
by imposing 



[(j){x,t),TT{y,t)\ ^ihS{x-y). 
Integration of y over the superfluid yields 

[kx),N] =-i, 
which in turn implies the uncertainty relation 



(73) 



A(f){x)AN > 



1 



If we assume that the number of particles in the super- 
fluid takes on a macroscopic value of order N w lO'^", a 
A7V of the order of ^/N implies a relative uncertainty of 
order 10"^" in the particle number and the phase. These 
numbers are comparable to the position and momentum 
uncertainties of a macroscopic object. The description of 
a macroscopic superfluid in terms of a classical field the- 
ory is therefore as appropriate as the classical description 
of any other macroscopic object. This is of course not in 



contradiction with the fact that Planck's constant H ap- 
pears in this effective field theory. We will see below 
that it manifests itself in a family of "quantum effects" , 
which are related to the compactness of the U(l) field 
(j){x). These effects require either a non-trivial topology 
or more than one superfluid, and are very similar for neu- 
tral and for charged superfluids. 



V. PHENOMENOLOGY AND THE HIGGS 
MECHANISM 

To begin with, however, let us consider the superfluid 
flow in a simply connected superfluid. The phenomenol- 
ogy depends strikingly on whether the fluid is charged or 
not. For a neutral superfluid, e* = 0, and the gauge fleld 
decouples completely. Even for a flxed set of boundary 
conditions, we have an infinite set of solutions for the 
superfluid flow, corresponding via 

m*Vs{x) — hV(j){x) 

to all possible choices of the phase field (f){x). In a simply 
connected superfluid, the flow will be vortex-free, i.e., 

V X Vs{x) = 0, 

and subject to boundary constraints, but apart from this, 
it only has to satisfy the continuity equation as an equa- 
tion of motion. 

The simplest example of a multiply connected super- 
fluid is a superfluid with a line defect, or vortex, along 
which the magnitude |4'(a;)| of the superfluid order pa- 
rameter vanishes. The phase (j){x) still has to be single 
valued everywhere in the fluid, but being a phase, its 
value may change by a multiple of 27r as we circumvent 
the line defect along a closed curve dS: 



as 



V4){x)dl = 27rn 



(74) 



where n is an integer. The angular momentum of each su- 
perfluid "particle" around the vortex is hence quantized 
in units of h. With Stokes theorem and the deflnition 



(jj{x) = V X Vs{x), 



(75) 



we may express this alternatively as quantization condi- 
tion for the vorticity 



/ u}{x)-nda = 
Js 



27rS7 

2m* 



(76) 



where n is a unit vector normal to the surface and the 
area integral extends over any open surface S which is 
pierced by the vortex once. The quantization of vortices 
in a superfluid is the simplest of the "quantum effects" 
alluded to above, where Planck's constant h enters in 
the phenomenology through the compactness of the fleld 
4>{x). (If (l){x) was not compact, we could eliminate h 
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completely from the effective theory by rescaling (f>{x) 

H{x).) 

Let us now turn to the phenomenology of a simply 
connected charged superfluid or superconductor, which 
displays the Higgs mechanism. The essence of the mech- 
anism is that the phase field 4){x) looses its independent 
significance in the presence of the gauge field. There are 
two ways of seeing this. The first is on the level of the 
equations of motion. We can simply choose a gauge such 
that (j){x) — everywhere in the fluid; any other choice 
of gauge can be brought into this gauge via (|77|l with 

he 

A(x) — (/-(x). 

The second is on the level of the effective Lagrangian. 
We may introduce a new vector field 



and p and J are given by H64() and H65|l . respectively. In 
principle, we could also obtain the continuity equation 
as the Euler-Lagrange equation for 0, but since L 
depends on derivatives of only through F^n and (j) is 
minimally coupled to A^, (|42|l is automatically satisfied 
by any solution of H79|l and H8()(l . This is consistent with 
the fact that has lost its independent significance due 
to the Higgs mechanism. 

For convenience, we choose the gauge 4>{x) = 0. Then 
and (ISni imply 



47re*(/9 - Us) 



Ane* Hs 



2ra*c^ 



m 



2m*c^ 



(81) 



(82) 



-A' = D., 



hd,,. 



-A, 



(77) 



In terms of this field, the superfluid Lagrange density 
(|63|l looks the same except that all terms containing (f) 
have disappeared. In particular, the terms quadratic in 
the derivatives of 6 in Cs have turned into a mass term 



2m* 



(78) 



for the vector field. The Maxwell Lagrange density and 
the Lagrange density for the uniform neutralizing back- 
ground charge take the same form with F^j^i, and A^ re- 
placed by F'^^ and A^, respectively, except for a total 
derivative or boundary term we discard. Thus the mass- 
less gauge field A^ is replaced by a massive vector field 
A^, while the Goldstone boson (f) disappeared. The to- 
tal number of degrees of freedom, however, is preserved: 
before, the massless vector field has two (the two helic- 
ity states of the photon) and the Goldstone boson one, 
while the massive vector field after the change of vari- 
ables has three degrees of freedom. In Sidney Coleman's 
words, "the vector field has eaten the Goldstone bosons 
and grown heavy" Q]. We will return to this issue after 
studying the phenomenology of the superconductor using 
the equations of motion. 

The Euler-Lagrange equation for , 



SL 



SL 



Sid^A,) SA, 



0, 



yields Maxwell's electrodynamics with electric charge 
density —e*{p~ rig) and current — e*J, 



V • £; -47re*(p- ris), (79) 

(80) 



1 Aire* 
V X B - -dtE^~^J, 
c c 



where 



Let us now restrict our attention to quasistatic phenom- 
ena, where we can neglect the time derivative terms. The 
analysis given below implies that this assumption holds 
for frequencies significantly smaller than c/Al, where 



47re*'^ns 



(83) 



is the London penetration depth. Then lf7^ - (|S^ reduce 
to 



t;2 A 2 V 



2m* 



V^A-V(VA) = -T^A\ 1- 



2m* 



(84) 



(85) 



Let us first look at the linear terms in these equations, 
i.e., the solution for infinitesimal and A. Under qua- 
sistatic conditions, (|80|l implies VJ — and with 182|) 
for infinitesimal fields VA — 0. The equations reduce to 



r-^ 1 
= - 

f 2 A/ 



(86) 



(87) 



i.e., we have electric screening in addition to magnetic 
screening, but with a screening length reduced by a fac- 
tor v/c. This leads us to conjecture that the dominant 
energy is the Coulomb interaction, which effects charge 
neutrality or p{x) ~ Us- We now simply assume that this 
is a valid approximation, and justify it a posteriori. Then 
(|HU implies 



2m* 



0, 



and reduces to 
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Taking the curl of this equation, we obtain London's 
equation jH ^ 



V X J = -B. 



(90) 



Under quasistatic conditions, we have again V J — and 
with VA = 0, which imphes that (|55|l reduces to 
(|S7|l . The solution of (|S7|l describes exponential screen- 
ing with penetration depth Al- If we have, for example, 
a superconductor which occupies the half-space x > 
subject to an external magnetic field B = Boy at the 
boundary x — 0, we obtain 



(91) 



where 



^0 — Al-Bq, "^0 



ATTm' 



■Bo. 



The screening of the magnetic field is known as the Meiss- 
ner effect. According to (I88f) . the vector potential implies 
an electrostatic potential 



B, 







-2x/\-L 



87re*ns 



(92) 



This potential allows us to verify the validity of our ap- 
proximation p{x) « Us- Substituting H92|l into H86(l . we 
find that the ratio of the neglected term to the terms kept 
is 



4^2 



(93) 



ie., the approximation is excellent. 

The electrostatic potential (|92ll is called the London or 
Bernoulli Hall effect p3.l23| . To understand its physical 
origin, it is best to rewrite H88() with H56|l for (l>{x) — in 
terms of the superfluid velocity: 



*^ 1 * 
-e $ -) — m Vg 

2 



0. 



(94) 



The electrostatic potential simply compensates the ki- 
netic energy contribution to the chemical potential, as 
required by H66|l with (f>{x) — 0. For stationary flow, 
this condition reduces to the requirement that the elec- 
trochemical potential /J.ci.chom. is constant across the su- 
perconductor. In practice, the London Hall effect can 
only be measured with capacitive contacts, as ohmic con- 
tacts are sensitive to the electrochemical rather than the 
electrostatic potential 24]. The effect furnishes us with 
an independent meaning of the superfluid density Ug or 
the effective mass to* , while under quasistatic conditions 
all other effects depend only on the superfluid stiff- 
ness rig/m*. The underlying theoretical reason is that 
the London Hall effect is a consequence of the correc- 
tions incorporated in the effective Lagrange density (|63|l . 
In these terms, the parameter m* enters by itself, while 



(apart from a total time derivative term irrelevant to the 
equations of motion) only the combination Us/m* entered 
in the previous approximative Lagrangian l|33|l with l|3U|l . 
^ and 

Since we have given a precise definition of the super- 
fluid velocity, it is legitimate to ask whether the London 
Hall effect (I94f) balances the Lorentz force, or, if not, 
what other forces balance it. The total electromagnetic 
force on a given particle with charge — e* in the fluid is 
given by 



(E+'^xB) 

^-dtA^S/^ + — X (V X A) 
c c 



(95) 



With (PH) and we obtain 



■ (dtVg + ^\/{vg)^ - X (V X Vg) 



m*{dtVg + {vgV)vg 



, dvg 
dt 



Thus the gradient V<I> of the electrostatic potential H94|) 
does not balance the Lorentz force -Vg x B, but both 
terms together account for the difference between local 
and substantial acceleration H70|l in the superfluid. This 
difference is significant when, for example, the flow is 
stationary but does not follow a straight line. 

Let us summarize how the Higgs mechanism manifests 
itself in the equations of motion. For a neutral superfluid 
with e* = 0, there are many solutions to 



J (x) = pVg 



_P_ 

TO* 



for fixed boundary conditions, corresponding to all pos- 
sible configurations for the phase field (t){x). These so- 
lutions reflect the existence of the massless sound mode 
described by the field (t>{x). For the charged superfluid, 
however, all the different configurations for (j){x) merely 
correspond to different choices of gauge; as far as the 
current or magnetic field distributions are concerned, all 
these solutions are always equivalent to one for which we 
have (f){x) = 0. Thus the field (j) does no longer describe 
an excitation. For a simply connected superconductor, 
it is only meaningful in that it assures gauge invariance, 
both on the level of the Lagrange density and on the level 
of the equations of motion. The solution of these equa- 
tions is physically (i.e., apart from the freedom to choose 
the gauge) unique for a given set of boundary conditions. 

We now return to the manifestation of the Higgs mech- 
anism on the level of the Lagrangian. For a simply con- 
nected superconductor, we have already seen that we may 
eliminate the phase field 4> if we introduce a new vector 
field A'^ according to lf77|l . The mass term (|75|) we find 
for A'^ may appear to violate gauge invariance, as mass 
terms generally do, and may, at first sight, to be taken 
as a signature of a spontaneously broken gauge invari- 
ance. For one thing, however, gauge invariance is not 
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violated. From the definition 17711 it is clear that the new 
field simply transforms as 

under a gauge transformation H27|l . The Lagrangian 
hence remains manifestly gauge invariant, and has to re- 
main gauge invariant, as it is the same Lagrangian as be- 
fore expressed in terms of different fields. Furthermore, if 
a symmetry is spontaneously broken, it is never violated 
on the level of the Lagrangian or the Hamiltonian, but 
only on the level of the ground state. 

In the literature, one sometimes finds the statement 
that "the gauge field acquires a mass" due to the Higgs 
mechanism. This is not exactly to the point, as it sug- 
gests that the massive vector field A'^^ is still a gauge field, 
while we have just seen that it is gauge invariant. In the 
case of a superconductor, we even know how to interpret 
the individual components of A'^ physically. According 
to 1771, (EZl, and 

(96) 

The Higgs mechanism hence does not imply that "the 
electromagnetic gauge field acquires a mass" , but only 
that we can describe the superconductor in terms of 
gauge invariant fields, that is, in terms of the chemical 
potential fi{x) and the superfluid velocity Vs{x). If we 
do this, we also have to express the Maxwell Lagrange 
density (|34|) in terms of /i and Vg. With 

we obtain for the total Lagrange density 

aire J 

(97) 

The Euler-Lagrange equations we obtain from H97(l for 
/i and Vs are equivalent to (f7i^ - l|H^ . and yield exactly 
the same solution as above. Writing the Lagrangian in 
terms of /i and Vg does not yield any practical advantage, 
but clearly illustrates that gauge invariance has become 
irrelevant — it is not broken, but has simply left the stage. 
Since all the fields are gauge invariant, (|97|l does not even 
provide a framework to think of a spontaneous violation 
of a gauge invariance. 

These considerations apply to every field theory which 
displays the Higgs mechanism. In any such theory, the 
Lagrange density is invariant under a global physical 
symmetry for a matter field, and invariant under a lo- 
cal gauge symmetry, which affects both the matter field 
and the gauge field. The global symmetry is "physical" 
as we can classify the states of matter according to their 



transformation properties, while the gauge symmetry is 
"unphysical" as gauge transformations have no effect on 
the states of matter, but only on our description of these 
states. In our example of a superfluid, charged or neutral, 
the global symmetry transformation is 

4>{x) -> -f A, (98) 

where A is independent of spacetime. This symmetry is 
spontaneously violated, which means that there are many 
degenerate ground states which map into each other un- 
der H98|l . For a neutral superfluid, we obtain a massless 
mode according to Goldstone's theorem. The situation 
is more subtle for a superconductor, as the matter field 
is coupled to a gauge field and the Lagrange density is 
also invariant under the gauge transformation (|27|l . This 
"unphysical" symmetry, however, seems to contain the 
physical symmetry as the special case 

he 

A{x) = --A. (99) 
e* 

The formal equivalence of the transformation (|98|) and 
(|27|l with (|99l) is at the root of the widely established 
but incorrect interpretation of l|98|) as a gauge transfor- 
mation, and in particular of the spontaneous violation 
of (|98|) as a spontaneous violation of a gauge symme- 
try. (This is presumably the reason why particle physi- 
cists like Steven Weinberg speak of "spontaneously bro- 
ken gauge symmetries" interchangeably with "the Higgs 
mechanism".) The problem here is that the equivalence 
is only formal. The gauge transformation (|99|l represents 
a transformation of our description, similar to a rotation 
of a coordinate system we use to describe a physical state, 
while the transformation H98|l corresponds to a transfor- 
mation of our physical state, like a rotation of a physical 
system. Clearly, a (counterclockwise) rotation of the co- 
ordinate system has the same effect on our equations as 
a (clockwise) rotation of the physical system we describe 
with these equations, but the transformations are all but 
equivalent. It is hence incorrect to refer to the spon- 
taneous violation of H98|l as a spontaneous violation of 
gauge symmetry. A gauge symmetry cannot be sponta- 
neously violated as a matter of principle. 

The difference between the "physical" symmetry (|98|) 
and the gauge symmetry H27|l can also be appreciated at 
the level of conservation laws. The former yields particle 
number (or charge) as a conserved quantity, according to 
(I42|l . while there is no conservation law associated with 
the latter. In the literature, (|98|l is often referred to as 
a global gauge transformation, and the conservation of 
charge attributed to gauge invariance. This view, how- 
ever, is not consistent. If one speaks of a global gauge 
symmetry, this symmetry has to be a proper subgroup 
of the local gauge symmetry group. The alleged global 
gauge symmetry hence cannot be a "physical" symme- 
try while the local gauge symmetry is an invariance of 
description, or be spontaneously violated while the lo- 
cal symmetry is fully intact. The difference between the 
global phase rotation H98|) and a global gauge rotation 
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(|99|l is even more at evident at the level of quantum 
states. The BCS ground state ^ is, for example, not 
invariant under H98|l. while it is fully gauge invariant, as 
we have seen in Section Ull 

The conclusions regarding the physical significance (or 
maybe better insignificance) of gauge transformations we 
reached here for superconductors hold for any field theory 
which displays the Higgs mechanism. 



VI. QUANTUM EFFECTS 

This discussion of the Higgs mechanism applies only 
to simply connected superconductors. If we have a non- 
trivial topology or more than one superfluid, the phase 
field 4> reassumes physical significance through its com- 
pactness, that is, the fact that its value is only defined 
modulo 27T. In these situations, we are not allowed to 
set (j){x) = in the equations of motion or eliminate it 
from the effective Lagrange density via H77|) or H96|l , as we 
would loose the information regarding the compactness. 
Since the phase field is multiplied by Planck's constant 
h whenever it enters in the Lagrange density, any effect 
due to the compactness of (j) will depend on h, and only 
exist for 0. Therefore we refer to them as "quantum 
effects". 

The simplest of these effects in superconductors is the 
quantization of magnetic flux, which is analogous to the 
quantization of vorticity in neutral superfluids. The ef- 
fect was predicted by London in a footnote in his first 
book |22.l2q almost a decade before BCS proposed their 
microscopic theory. Consider a macroscopic supercon- 
ductor with a hole in it, which may either be a hole in 
the superconducting material or a line defect or vortex in 
the superconducting order parameter. Like in the case of 
a vortex in a neutral superfluid, the phase field 4) has to 
be single valued everywhere in the superconductor, but 
may change by a multiple of 2tt as we circumvent the 
hole or defect along a closed curve dS: 



V(i){x)dl = 27rn, 



(100) 



dS 



where n is an integer. We now take dS well inside the 
superconductor, that is, separated at each point by a 
distance much larger than the penetration depth Al from 
the hole or defect. Then, according to the Meissner effect 
or our derivation of London's equation above, which still 
applies locally, the superfluid velocity 

1 / e* \ 

= — A] 

m* \ c I 

has to vanish along 95, and (|100|l implies 

I A(x)dl^ [ B{x)-nda= — -n, (101) 
JdS Js e* 

where we have used Stokes theorem once more. The mag- 
netic flux through the hole or vortex is hence quantized in 



0— ' 

vortex 



(a) 



(b) 




FIG. 2: Phase slippage: A vortex moving in a superfluid in- 
duces a transverse gradient in the (electro) chemical potential 
by dragging a branch cut in the phase of the order parameter 
with it. 



units of hc/e*, which for e* — 2e is half of the Dirac flux 
quantum. Note that the vorticity H75|l is not quantized 
in a superconductor. 

We now review two further quantum effects, which 
are similar in neutral and charged superfluids; as in our 
derivation of the effective theory above, the equations 
for the latter case contain the former as the special case 
e* — 0. One of the effects is phase slippage 0- Consider 
two points 1 and 2 in a superfluid, which are connected 
by a vertical path (see Fig.l^^). Now imagine we adia- 
batically move a vortex from the very far left across the 
path to the very far right. This process yields a difference 
in the electrochemical potential between the two points, 
which is according to (|68|l given by 



'dt (^hW<j>{x) + ^Mx)) dl. (102) 

where the line integral is taken along the path between 
the points. The time integrated difference in the elec- 
trochemical potential is hence given by the difference be- 
tween the line integral at the end of the process and the 
line integral at the beginning. Let us first consider the 
case of a neutral superfiuid, i.e., e* = 0. The line integral 
of V(f> will have changed by 27r, as the difference in the 
paths is topologically equivalent to encircling the vortex 
once (see Fig.lJb). Alternatively, we may say the vortex 
has dragged a branch cut across the path. (We assume 
that at the beginning and the end, the vortex is so far 
away from points 1 and 2 that we can neglect its influence 
on the line integral.) If we now have a continuous flow 
of vortices across the path, the line integral will pick up 
a contribution of 27rfi, from each of them, and we obtain 
a chemical potential difference 



(103) 



where (9tiVv)av is the average rate of vortices crossing the 
path. 
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Let us now turn to the case of a superconductor, where 
we assume that during the entire process the distance be- 
tween the vortex and either of the points 1 and 2 is much 
larger than the penetration depth. Since the hnc inte- 
gral we obtain when encircling a superconducting vortex 
along a circle well inside the superconductor is zero, 



-R»Ai 



hV4>{x) 



A{x)\dl = 0, 



(104) 



we do not obtain a difference in the electrochemical po- 
tential as we move an isolated vortex carrying a magnetic 
flux quantum across the path. So, at first sight, it may 
appear as there is no phase slippage effect in supercon- 
ductors. The situation just described, however, is not the 
general one, as we dragged a unit of magnetic flux with 
the vortex from the very far left to the very far right. 
This produced a Hall effect which exactly canceled the 
phase slippage effect. If we consider a situation where 
we have a large, almost uniform magnetic field and an 
Abrikosov vortex lattice or liquid in which the distance 
between the vortices is much smaller than the penetra- 
tion depth, and we have a flow of vortices across the path, 
the magnetic field will remain to a reasonable approxi- 
mation unaffected by the flow and we recover IjlQ^II for 
the electrochemical potential difference. The voltage we 
measure between the two points is then given by 



U 



2e 



(105) 



This voltage is known as the Nernst effect in supercon- 
ductors. 

The last and possibly most striking quantum effect we 
review is the Josephson effect |23|. Consider two super- 
fluids or superconductors Si and S2, which are weakly 
coupled, say by a narrow constriction for superfluid par- 
ticles or a tunneling barrier for Cooper pairs. The only 
requirement for the effect is that there is an energy asso- 
ciated with the weak link, which depends on the (gauge 
invariant) phase difference A0 between two points 2 and 
1 in superfluids S2 and Si: 



'junction 



with 



A<t> = <t>{2)-4>{l) + ^^j^ A{x)dl, 



(106) 



(107) 



where the line integral is taken along the path the super- 
fluid particles take when they move from one superfluid 
to the other. Note that A(/) divided by the distance be- 
tween the points 2 and 1 is just the discrete version of 
the gradient term 

hV(t){x) + —A{x) 



given by its square. We assume that -Ejunction is likewise 
minimal for A(j) = 0, which implies that the first term in 
a Taylor expansion around this minimum is quadratic in 
A(f). In the case of the junction, however, this term is not 
sufficient. Since (f> is only defined modulo 2tt, /(A0) has 
to be a periodic function of A(j). Josephson has shown 
that to a reasonable approximation, it is given by 



/(A(/.) = -^ocos(A(/.). 



(108) 



Let us now assume a situation where both macroscopic 
superfluids are in a state of equilibrium, but the phases 
are not necessarily aligned relative to each other. Then 
only the energy stored in the junction depends on the 
phases of the superfluids, and the "characteristic equa- 
tion" (|7^ becomes for superfluid S2 



hdtN2 = 



'junction 



(A0) 



dm 



(109) 



where A'2 is the number of superfluid particles or Cooper 
pairs in S2 . (We would also obtain a similar equation for 
Si, but since we assume Ni + N2 — const, and injunction 
only depends on (/i(2) — 0(1), it does not provide any addi- 
tional information.) The particle current from superfluid 
Si to S2 is hence given by 

1 disjunction (A0) 1^ ■ / \ j.\ /11n^ 

Jl^2 = j:^^^ = j:EoMm (110) 

On the other hand, since the other "characteristic equa- 
tion" (|66|l holds for each superfluid, 

hdt{m - m) = -{m - + e*{^2) - $(1)). 

(Ill) 

If we add 



dt 



A{x)dl 



to both sides of Hlll(l . wc obtain 

HdtA^ = - nil)) -e* J E{x)dl - -Aa..,.,_, 

(112) 

or, if we take A/i„,^i-,„,-„ time independent, 



Substitution into (|110() yields 
Eo 



where 



Ji^2 = sin(27ri'i -|- A(/)o) 



A^ol.cho 



(113) 



(114) 



we already encountered in the Ginzburg-Landau free 
energy, where the magnetic energy (|25f) was essentially 



is the Josephson frequency. This implies that if the elec- 
trochemical potential is equal for both superfluids, we 



find a DC particle current depending on the initial align- 
ment of the phases. If there is a difference in the poten- 
tial, however, the current will oscillate with frequency v. 
This is called the AC Josephson effect. The effect exists 
for both neutral and charged supcrfluids, but it is much 
easier to measure in a superconductor, as we can realize 
a difference in the electrochemical potential by applying 
a voltage U across the junction, AyLt„,^i,cm. = — 2eC/, and 
easily measure oscillations in the electrical current. 

Note that the Josephson effect, so astonishing its phe- 
nomenology may be, follows through the "characteris- 
tic" equations of superfluidity directly from the fact that 
there is a broken symmetry in superfluids and that the 
compact phase field which labels the different degenerate 
ground states is the field conjugate to the density in the 
superfluid. The other assumption we made in this arti- 
cle, the assumption that both current and momentum are 
carried by the same species of particles in a superfluid, 
was not required to explain any of the quantum effects. 

APPENDIX 

In this appendix, we derive the position space wave 
function Q by projecting the BCS state onto a flxed 
number of pairs N . The (unnormalized) BCS state may 
be written 
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k 

k 

= exp(e*'^&^) |0). 
The pair creation operator 6^ is given by 

where (p{x) is given by ((HJ. If we now project out a state 
with N pairs according to we obtain 

\M= I #e-^^*exp(e^^&t)|o) = |.(6t)^|o), 

which is (up to a normalization) equivalent to 
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